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Abstract. A scheme for fast, compact, and controllable acceleration of heavy 
particles in vacuum is proposed, in which two counterpropagating lasers with variable 
frequencies drive a beat-wave structure with variable phase velocity, thus allowing 
for trapping and acceleration of heavy particles, such as ions or muons. Fine 
control over the energy distribution and the total charge of the beam is obtained via 
tuning of the frequency variation. The acceleration scheme is described with a one- 
dimensional theory, providing the general conditions for trapping and scaling laws for 
the relevant features of the particle beam. Two-dimensional, electromagnetic particle- 
in-cell simulations confirm the validity and the robustness of the physical mechanism. 
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1. Introduction 

Compact acceleration of charged particles has become a central research topic due to 
the recent developments in the technology of ultraintense lasers in the infrared (IR) 
optical regime. So far, the main research directions in the field have been electron 
acceleration in laser-driven plasma waves PQ [5], and proton or light-ion acceleration in 
laser-solid interactions (cf. [6] [10] and, for a review, [TH [12]), which is of particular 
interest for possible applications to the treatment of tumors with hadron beams |13j-|17j. 
These acceleration schemes are indirect, in the sense that the acceleration is provided 
by laser-driven field structures in plasmas. In the case of electrons, the availability 
of relativistic radiation intensities (/ > 10 18 W/cm 2 with IR lasers) also motivated a 
number of proposals for direct, all-optical acceleration in vacuum by the electromagnetic 
(EM) field of one or more laser pulses |18]-[30j: in the case of ions, for which relativistic 
radiation intensities (/ > 10 24 W/cm 2 with IR lasers) are still beyond the limits of 
current technology, direct acceleration has been mostly unexplored (only very recently, 
the possibility of direct ion acceleration by tightly focused, radially polarized lasers has 
been investigated [3T]). 

In this paper, a new concept for direct, all-optical acceleration in vacuum is 
presented, which allows for fast, compact acceleration of heavy particles, including 
protons, heavy ions, and muons. The technique exploits two counterpropagating EM 
waves with variable frequency in order to drive a beat-wave structure with controllable 
phase velocity. If the frequency variation is chosen appropriately, the ponderomotive 
potential associated with the beat wave provides longitudinal trapping and long-term 
particle acceleration. The main characteristics of the method are the following: (i) it 
can work with "modest", nonrelativistic radiation intensities (e.g., / <C 10 24 W/cm 2 
with IR lasers), provided that the available electromagnetic energy is sufficiently high 
(typically, > 100 J); (ii) it does not involve any complex laser-matter coupling; (iii) it 
exploits an extremely simple and robust physical mechanism, operating in a test-particle 
regime, as a conventional accelerator, offering the possibility of direct, efficient control 
over the relevant beam features (mean energy, energy spread, and total accelerated 
charge) via regulation of the fundamental laser parameters (intensity, duration, and 
frequency variation); (iv) it can be adapted to any source of charged particles (including 
beams, and tenuous gases and plasmas), serving either as a pre-accelerator, creating 
an energetic beam from heavy particles with low energy, or as an energy booster, 
accelerating particles from a beam with a given energy. 

The main advantages of the acceleration technique, controllability, robustness, and 
simplicity of the physical scheme, are obtained through two technological requirements 
[32] : (i) the method typically requires long, energetic laser pulses to maintain the 
accelerating structure for sufficiently long times; (ii) the method typically requires 
significant frequency excursions to achieve high energy gains. If these technological 
requirements can be met, the technique can represent an alternative option whenever 
beam quality (e.g., a low energy spread) and a high degree of control on the particle beam 
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are critical (the pertinence of the present technique with respect to other laser-based 
techniques depends on the specific application). Furthermore, the present acceleration 
scheme could also be an option in applications for which no other laser-based acceleration 
techniques are currently available [33] (e.g., ultrafast acceleration of muons during, or 
immediately after, the muon cooling process in a muon collider or a neutrino factory 
machine [31] [IT]). 

The Paper is organized as follows: in Section [21 the basic physical mechanism of 
the acceleration scheme is described with a one- dimensional theory; in Section [3j the 
resonant solutions of the system and the conditions for particle trapping are analyzed; 
in Section HI basic scaling laws for the properties of the particle beam are provided; 
in Section [5j numerical results from self-consistent simulations in two dimensions are 
presented. Finally, in Section El the conclusions are stated. 



2. Physical mechanism 

The essential aspects of the method are illustrated with a one- dimensional, single- 
particle theory, in which the amplitude of both lasers is constant, and only one of the 
beams has variable frequency. For the sake of simplicity, the Lorentz reference frame 
where the initial velocity of the particle is zero, and the laser frequencies at x = are 
initially the same, is adopted. The first laser beam, propagating from left to right in 
the x direction, is described by the vector potential Ai = Ai{cos(9) sin[fco£ + <t>{£)]&y + 
sin(6 l ) cos[&o£ + <KO]^z}? where 9 determines the polarization type (e.g., 9 = n/2 for 
linear polarization along z, 9 = 7r/4 for circular polarization), £ = x — ct, ko = 2tt/\o 
is the wavenumber at £ = 0, and <j) is an arbitrary function, such that its derivative 
vanishes for £ = 0. The second beam propagates in the opposite direction, with vector 
potential A 2 = A 2 {cos(6 l ) sin[A;o(— x — ct)]e y + sin(#) cos[/co(— x — ct)]e z }. By resorting 
to the conservation of the transverse canonical momentum, the x component of the 
equation of motion for the particle is written as 

where q and M are the particle charge and mass, respectively, and 7 is the Lorentz 
factor. If Aj = <C 1, for j = 1,2, the longitudinal displacement due to the quiver 
motion can be neglected (here, Aj = ^CLj, where m and e are the electron mass and the 
elementary charge, and aj = ^—^ is the usual normalized vector potential, corresponding 
to the peak transverse momentum of an electron in the laser field). Hence, as long as 
the ratio between the low and high beat-wave frequencies is negligible along the particle 
trajectory x p (t) - i.e., as long as |ti> p /fl p | 1, where u p — (c — f p )y — k v p and 
Q p = (c — v p )y + koc, with v p = - the equation of motion can be averaged over the 
fast time scale to yield, in dimensionless form (t = k ct, x = k x, and p x = j^), 

dp x A 1 A 2 d 



di 27 dx 



cos [2x + <j>{x - t)] , (2) 
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where 7 = {1 + p 2 x + A 2 /2 + A\/2 + A x A 2 cos[2x + (f)(x - i)}} 1 ' 2 now denotes the 
average Lorentz factor [22] ■ Equation ([2]) can be obtained from the Hamiltonian 
J4?(x,p x ,t) = {l+p 2 x + A 2 1 /2 + A 2 /2 + A l A 2 cos[2x + (j)(x-t)}} 1 / 2 . Hence, the variation 
of the particle energy, = is given by 

¥ = ^I COS[2£ + 0(£ - £)] - (3) 
If there is no frequency variation (<// = 0), J$? is time independent and 7 remain constant 
(in agreement with the Lawson- Woodward theorem [HI [43], Hi]). In the presence of 
frequency variations, the phase velocity of the ponderomotive beat wave varies in time, 
allowing for wave-particle energy transfer. If <p is chosen appropriately, the particle 
is trapped by the accelerating beat wave: since, for a trapped particle, the condition 
|c<jp/fip| <C 1 is always satisfied and Eq. ([2]) is valid for every t, the particle can be 
accelerated continuously and over long times. 

3. Phase-locking and particle trapping 

For a particular choice of the function 0, equation (TSJ) admits a resonant solution, X(t), 
such that 2X(t) + 4>[X(t) —i] = 4>o- In this solution, the particle is exactly phase-locked 
to the beat wave, i.e., its trajectory is identical to the trajectory of the point of the beat 
wave having constant phase 0o> here denoted with x^i). The resonance occurs when 
depends on £ = x — t as 

0(O = 0o-f+^log(l + 2/if), (4) 

where /1 = AiA 2 sin(0 o )/7o x , with 7q ± = 1 + A\/2 + A 2 ,/2 + A X A 2 cos(0 o )- The resonant 
trajectory can be expressed in parametric form, as a function of the proper time r 
(defined as dr/di = 7 _1 ), as 

X(t) = ^-[fi T^ r/2 - 1) - log(l - ii q t)}, (5a) 

t{r) =^[/vr(l-/W2)-log(l-// T)], (5b) 

with /i = 7o x /i. Accordingly, the particle energy increases as j(t) = 7o x [l + 
(1 — fi r) 2 ]/[2(l — (jlot)]. For t yU -1 , the energy grows quadratically in time as 
7o x (l + ^ 2 i 2 /2). For i > fi~\ the energy scaling depends on the sign of /z: if 
/i > (the particle and the variable- frequency laser are copropagating), the frequency 
must be increased to maintain phase-locking and the energy grows exponentially as 
7o x /2 exp(2/xt — 1/2); if [i < 0, the laser frequency must be decreased and the energy 
scales with the power law 7o x (|/i|t) 1 / 2 . 

In order to determine the condition under which the acceleration mechanism is 
effective, the stability of the particle motion has been investigated. The analysis reveals 
that stable solutions, and, consequently, the occurrence of particle trapping, exist 
whenever the frequency variation is sufficiently slow, independently of the particular 
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shape of 0. Expressing Eq. ([2D in terms of the phase difference between the particle 
and the wave, ip = 2[x — x^, (i)], one obtains 

d 2 ^ d 



dr 2 dtp 



uy>,r), (6) 



where 



U(^,t) = 2i 1 i 2 cos{^ + 2^Jt(r)] +0[4(r) +^/2]} +2a 0o (r)^, (7) 

with ^ (r) = x^ [t(r)] ~ K T ) ancl «</>o( r ) = ^^o[*( r )]- Tne quantity U plays 
the role of an effective potential and its shape determines the possible occurrence of 
particle trapping, which is allowed only if the frequency variation is smooth enough to 
guarantee that U presents local minima. For trapped particles, \if)\ <C 2^ and U can 
be approximated as 2A\A 2 cos{[l + (j)'{^ ) /2\ip + O } + 2a^ ip, leading to the necessary 
condition for particle trapping 

|a^| < a M = A X A 2 [l + 0' (l^ /2] , (8) 

where a M is the maximum value of the ponderomotive force. The physical interpretation 
of Equation (JSJ) is that trapping is allowed only if the maximum ponderomotive force 
of the beat wave is greater than the inertial force associated with the beat-wave 
acceleration. The limit in which \ct^ \ = a M corresponds to the resonant solution of 
Equations fl5aj) and (j56| . in the particular case O = ±tt/2. When is given by (TJJ 
and 1 0o I 7^ 7r/2, regions where trapping is possible always exist, although the solution 
is stable for cos(0 o ) < and unstable for cos(0 o ) > 0, with -0 = corresponding to 
the bottom and the top of the potential well, respectively. The initial relative width of 
the potential well, ~ 1 — - arcsin[a^ o (0)/a M ], provides an estimate for the trapping 
efficiency. As an example, the typical trajectories of three test particles are shown in 
Fig. [U along with the evolution of the beat-wave ponderomotive force. 

The generality of the trapping condition guarantees that the resonance criterion 
fll]) can be relaxed, provided that is sufficiently smooth, opening the possibility of 
controlling the final particle energy and the total accelerated charge by tuning the 
frequency variation. An important choice for 0, particularly relevant in experiments, is 
0(0 = 0o + o~i 2 , corresponding to a linear chirp. In this case, the beat- wave trajectory, 
x^ (t) =£-[l-(l-2at) 1 / 2 ]/(T, is such that a^(i) = -a{l ~2ai) 1 ' 2 /\l - 2(1 - 2ai) l l 2 } 2 . 
For t — 0, trapping regions exist if |er| < A\A 2 . The energy of the trapped particles 
grows approximately as [2(1 — 2ai)~ 1 / 2 — (1 — 2ot)~ 1 ]~ 1 / 2 . Depending on the sign of a, 
the trapping region evolves in time differently: if a < (trapped particles move from left 
to right), it becomes wider (cf., movie movl in [15]); if o > (trapped particles move 
from right to left), it becomes narrower, until disappearing. Thus, long-term stable 
acceleration can be achieved only when the accelerating particles and the chirped laser 
are copropagating. 
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Figure 1. Space-time evolution of the dimensionless ponderomotive force — J| cos[2x+ 
4>(x — t)] (color scale, red and blue denoting accelerating and decelerating regions, 
respectively), and trajectories of a resonant particle (solid line), a non-resonant, 
trapped particle (green markers), and a non-resonant, untrapped particle (yellow 
markers), for a case where (f> is given by Equation ((4]), A\ A 2 — 5 x 10~ 4 , and 0o = 57r/6. 
The inset represents the effective potential U at t = 0, gray marking the region where 
trapping is possible, and the corresponding initial positions of the three particles. 



4. Estimates for the basic beam properties 

In the nonrelativistic regime, the maximum energy gain, Ae M , can be estimated 
as Ae M [MeV] « 0.8 q 4 [e] M- 3 [amu] ^[lO^W/cm 2 ] J 2 [10 20 W/cm 2 ] A 2 [/mi] AT 2 [ps]; 
the corresponding acceleration distance scales as Ax [/mi] « 6 q 2 [e] M~ 2 [amu] 
^[lO^W/cm 2 ] J 2 1/2 [10 20 W/cm 2 ] A [/H AT 2 [ps]. By varying the chirp rate a, the 
energy gain Ae can be controlled, according to Ae = (o-/<j M ) 2 Ae M . Since the trapping 
efficiency decreases when increasing \a\, the chirp rate must be chosen in order to obtain 
the desired trade-off between energy gain and total accelerated charge. 

By analyzing the range of validity of the model, conclusions can be drawn for more 
general situations, in which spatial distributions of charged particles are considered 
and the effects due to the finite-size of the laser pulses are taken into account. In the 
presence of a spatial charge distribution, the scheme remains effective if the particle 
density n is low enough that the space-charge field does not perturb the EM beat-wave 
structure. Assuming that a single particle species is present, this establishes the criterion 
n[10 19 cm- 3 ] < g- 1 [e]/ 1 1/2 [10 20 W/cm 2 ]/ 2 1/2 [10 20 W/cm 2 ]. When taking into account the 
finite size of the laser beams, both the longitudinal and the transverse dependence of the 
laser envelope must be considered. Since long pulses (e.g., with duration in the range 1 
ps - 1 ns) are required to obtain significant energy gain, A\ and A2 can be considered 
as slowly varying functions of t. Hence, intensities high enough for particle trapping 
exist only for a fraction of the overlap between the laser pulses, and the duration of the 
acceleration process is limited to a fraction of the pulse duration, determined by the 
specific chirp law. For similar reasons, the transverse cross section of the accelerating 
region is on the order of tiWq, where Wo is the overlapped spot size of the lasers, and 
the length of the acceleration process is always limited to distances on the order of the 
Rayleigh length, Z R = ttWq/Xq. Therefore, the total accelerated charge Q scales as 
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Q[pC] ~ 16^ tr g[e] n[10 19 cm~ 3 ] Wq [/xm]/Ao[/xm], where rj tT is the fraction of particles 
trapped by the beat wave. The finite transverse size of the lasers also causes a slow 
outward drift due to the transverse ponderomotive force of the lasers. If the intensities 
of the two lasers are comparable, the ratio between transverse and longitudinal momenta 
scales as ^ ^ Transverse focusing and trapping can be obtained using laser pulses 
with appropriate shape [16]. 

The above scalings indicate that the acceleration scheme described here does allow 
for compact, fast acceleration of heavy particles. In the case of protons, proof-of- 
principle experiments producing monoenergetic beams in the 100 keV - 1 MeV energy 
range could already be performed using kJ-level, suitably chirped pulses, with frequency 
excursions on the order of 5-10% of the central frequency (cf. [32J), which are potentially 
available with current laser technology. Furthermore, the technique could open the way 
to ultrafast acceleration of low-energy muons [33], for which multi-MeV energy gains, 
on a ps time scale, could be achieved using lasers with moderate energy (on the order 
of 1 J) [17]; such a possibility could be critical to minimize muon-decay losses. 

5. Numerical simulations 

The effectiveness of the method has been confirmed with EM particle-in-cell simulations, 
using the OSIRIS 2.0 framework [18], accounting self-consistently for space-charge and 
propagation effects. Here, a selection of results is presented from a two-dimensional 
(2D) simulation, in which the source of particles to be accelerated is provided simply by 
a slab of hydrogen, with thickness 75 /im and proton density n = 5 x 10 16 cm~ 3 . The 
computational domain is rectangular [L x = 300 /im (5 x 10 4 cells) and L y = 50 /xm (100 
cells)], with open-space boundary conditions and periodic boundary conditions in the x 
and y directions, respectively; the left face of the hydrogen slab is located at x = 40 /im. 
Two approximately gaussian laser pulses, with spot size Wq = 10 /im and duration 4.2 
ps, are focused on the mid plane of the slab. The pulse propagating from right to left has 
peak intensity J 2 = 8.5 x 10 20 W/cm 2 , fixed frequency, and central wavelength A = 800 
nm. The pulse propagating from left to right has peak intensity I± = 1.3 x 10 21 W/cm 2 
and its phase varies quadratically as k £ + cr(£ — £o) 2 ? where £ indicates the distance 
from the pulse center, a = —2 x 10~ 5 fc 2 is the chirp coefficient, and £o — 6.25 x 10 3 kg 1 
is the point of the pulse where k = ho, which need to be specified when using finite- 
length pulses. Field ionization of the neutral atoms [H] is calculated by using the 
Ammosov-Delone-Krainov rates [51)]. The laser pulses are launched from the opposite 
x boundaries of the computational domain and hit the left face of the hydrogen slab 
at the same time. After 1.8 ps, the trapping process begins and 4.2% of the protons 
initially within the focal region of the lasers are extracted from the bulk distribution 
and accelerated (cf., movie mov2 in Ref. [IS]). As the acceleration continues, these 
protons form a bunch that separates from the main distribution both in space and in 
momentum, as illustrated in Figs. [2] and [3J After 6 ps, the proton beam has a mean 
energy of 12 MeV with a 7.5% energy spread, approximately 60 /im long and 20 /im 
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Figure 2. (Color online) Proton distribution in the x-p x phase space, after 6 ps. The 
inset shows the corresponding kinetic-energy spectrum. 



wide, with a transverse momentum spread on the order of 10~ 2 Mc. Much lower energy 
spreads, down to 1% and below, can be obtained by employing longer pulses with the 
same energy and lower intensities. The expansion of the remaining plasma bulk, visible 
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Figure 3. (Color online) Mean kinetic energy of the protons as a function of position, 
after 6 ps. 



in Fig. [2] and in mov2 in [15] , is driven by the space charge that forms because of the 
rapid expansion of the hot electrons produced by the lasers (for the radiation intensities 
considered here, the electron dynamics in the beat wave is dominated by stochastic 
heating [511 [52]). These results confirm the potential of the proposed technique for 
direct acceleration as a method for producing high-quality beams of heavy particles in a 
controlled way. Other simulation results, to be presented in future works, suggest that, 
for higher plasma densities, the space-charge field can affect the dynamics and, under 
appropriate conditions, it can even improve both trapping and acceleration. 



6. Conclusions 



Trapping and acceleration of heavy particles in a purely optical structure, generated by 
two counterpropagating laser beams with variable frequency, has been demonstrated. 
Relying on an extremely simple and robust physical mechanism, the method proposed 
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here offers unique capabilities of detailed control over the main features of the particle 
beam. This makes the new acceleration technique a promising candidate for any 
application where the controllable production of high-quality beams of heavy particles in 
compact machines is critical, and for applications where ultrafast acceleration is crucial 
(e.g., the production of muon beams in a muon collider or a neutrino factory machine) 
Pi- 
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